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ABSTRACT 

G349.7+0.2 is a supernova remnant (SNR) expanding in a dense medium of molecular clouds and 
interacting with clumps of molecular material emitting gamma rays. We analyzed the gamma-ray data 
of Large Area Telescope on board Fermi Gamma Ray Space Telescope and detected G349.7+0.2 in the 
energy range of 0.2—300 GeV with a significance of ^13cr showing no extended morphology. Modeling 
of the gamma-ray spectrum revealed that the GeV gamma-ray emission dominantly originates from 
the decay of neutral pions, where the protons follow a broken power-law distribution with a spectral 
break at ^12 GeV. To search for features of radiative recombination continua in the eastern and 
western regions of the remnant, we analyzed the Suzaku data of G349.7+0.2 and found no evidence 
for overionized plasma. In this paper we discuss possible scenarios to explain the hadronic gamma-ray 
emission in G349.7+0.2 and the mixed morphology nature of this SNR. 

Keywords: ISM — ISM: clouds — ISM: individual object (G349.7+0.2) — ISM: supernova remnants 
— X-rays: ISM — gamma rays: ISM 


1. INTRODUCTION 

The Galactic supernova remnant (SNR) G349.7+0.2 
has been long suggested to be on the edge of the Milky 
Way ( ^22 kpc), following t he Hi absorption measure¬ 
ments (Caswell et al. 19751) and the OH (1720 MHz) 
maser observations (lb rail et al.||l996|). Recently, based 
on a better understanding of the Galactic Center’s struc¬ 
ture and using the 1420 MHz radio continuum and 21 cm 
Hi emission data sets from the Southern Galactic Plane 


Survey, Tian & Leahy (2014) estimated the distance of 
G349.7+0.2 to be ^11.5 kpc. 

In the radio waveband, the central spectral index value 
of —0.47 ± 0.06 suggests that G349.7+0.2 is a shell-type 
SNR. However, due to the enhancements towards the 
eastern and southern parts of its 2.'5 diameter shell, it was 
found to be an unusual shell-typ e SNR lacking a ce ntral 
cavity and a ring-like structure (Shaver et al. 1985). 

In X-rays, G349.7+0.2 w as first detected Tp the Galac- 
tic plane survey of ASCA ( Yamauchi et a l. 1998). The 
analysis of the ASCA data (j Slane et al.[ 


showed 

that the spectrum was well tit by a single-temperature 


non-equilibrium ionization (NEI) model and they esti- 
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mated the ag e of G349.7+0.2 to be ~2800 yr. |Lazendic 
et al. (2005) analyzed the Chandra data of the whole 
S1NR as well as six spectral regions of the remnant. The 
spectrum of the whole SNR was fit by two thermal com¬ 
ponents: a high-temperature plasma (~1.4 keV) in NEI 
with an enhanced Si abundance and a low-temperature 
plasma (^0.8 keV) in ionization equilibrium. Unlike the 
whole SNR spectrum, the spectra from individual SNR 
regions were well fit with a single-temperature thermal 
model. The Chandra observations also revealed a com¬ 
pact central object inside the SNR shell, called CXOU 
J171801.0—372617. The X-ray morphology was found to 
be very similar to the radio morphology, thereby giving a 
hint of an expansion into a region with large de nsity gra¬ 
dients as it is shown by the 12 CO obs ervations (Reynoso 
fc Magnum||2001 Dubner et al.||2004 ). 

Five OH maser spots were found in the velocity range 
of 14.3 and 16.9 km s -1 inside G349.7+0.2 showing clear 
evidence of interaction between the SNR and the moleo 
ular 

at the brightes t 

(Brogan et al.p OOO). The velocities of these OH maser 
spot s are m agreement with the velocity of the sho eked 


clouds (MCs) (|Frail et al.| 19961). The magnetic field 
he brightest maser is measured to be 350 ± 50 ptG 


MC ([Reynoso fe Magnum||20Qll |Du bner et al.|2004|) and 
the v elocity of the far 3 kpc arm (Dame fc Thaddeus 
2008). 


Molecular line transitions found at similar velocities 
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of those of the OH masers and the shock-excited near- 
infrared H 2 emission hints that the SNR-MC i nteractio n 


is happening toward the center of the SNR (Lazendic et 
2010 Tian & Leahy|2014). By looking at the CO 
aata and the re sults of the HI absorption analysis, |Tian 


& Leahy (2014) concluded that the clouds with a velocity 
of 16.5 km s _r are behind G349.7+0.2 and that they are 
shocked by the SNR, where the interaction gives rise to 
near-infrared H 2 emission and the OH maser emission. 

It has been suggested that all maser-emitting (ME) 
SNRs could be of the mixed-morphology (MM) class 
Rho & Petrel |1998| lY usef-Zadeh, Wardle, fc Robertsl 
by [Lazendic et al. (2005) that the 
X-ray morphology of G349.7+0.2 is not centrally peaked. 
It was suggested that this morphology can be explained 
by the projection effect model by Petruk (2001), which 
assumes a shell-like evolution of an SNR at the edge of 
an MC, where the ambient density gradient does not he 
in the projection plane. So it is highly likely that the 
G349.7+0.2 is expanding into a density gradient with an 
angle of ^45°, rather than expanding along the direction 
of the density gradient. 

In gamma rays, MM SNRs interacting with MCs were 
among the first SNRs detected by the Large Area Tele¬ 
scope (LAT) on board Fermi Gamma Ray Space Tele¬ 
scope (Fermi- LAT). The gamma-ray luminosities of MM 
SNRs were found to be much higher than that of other 
detected SNRs, i.e. ^10 35 —10 36 erg s -1 for IC443, W28, 
W51C, W44, W49B, and 3C391 |Abd~o et al.|[20loiTb 


2009a 2010c d Castro & 


Slane||20l0j|Ergin et al.||2014p. 

SNRs with IVICs were clearly 


Interactions of these MM 

shown by the detection of 1720 MHz OH masers ( Yusef- 
Zadeh, Uchida, fc Roberts 1995f [Frail et al.|| 19961 Ureen 


et al.||1997[ |Claussen et al.j 1997 


Hewitt & Yusef-Zadeh 


al. 

2007 Reach, Rho & Jarrett 2005 

. 

lor G349.7+0.2, Castro & Slane 

(2010| 


reported a 
the position 

of G349.7+0.2 from the data collected by Ferrm-LAT. 
No evidence of spatial extension was found and a power- 
law fit to the gamma-ray data gave a spectral index of 
T = —2.10 zb 0.11. The shape of the gamma-ray spec¬ 
trum showed a steepening of the spectral index after a 
few GeV. So, a fit to an exponential power-law function 
resulted in a slightly better fit with a cut-off energy of 
~16.5 GeV a nd spectral index o f T = —1.74 =b 0.37. In 
TeV energies, Abramowski et al. ([2014 ) reported the de¬ 
tection of G349.7+0.2 at 6.6cr and found the TeV flux at 
roughly 0.7% of the Crab Nebula’s gamma-ray flux with 
a spectral index of 2.8 dt 0.27 sta t ± 0.20 syst . 

MM SNRs interacting with MCs are primary tar gets 
for the dete ction of gamma rays of hadronic origin (Slane 
et al. 2014), where two gamma rays are produced from 
the decay of a neutral pion (7r°) created in a proton- 
proton interaction while the SNR shock passes through 
a dense molecular material. The gamma-ray spectra of 
these SNRs fall steeply below 250 MeV and at energies 
greater than 1 GeV they trace the p arent proton energy 
distribution (Ackermann et al.||2013). 

In addition to the fact that G349.7+0.2 is an ME 
SNR and shows a high gamma-ray lumin osity < 

f r ^ J 10 35 


erg s 1 ) in the GeV energy range (Castro fc Slane|2010) 


of G349.7+0.2, namely, the gamma-ray spectrum be¬ 
ing dominated by the hadronic emission component and 
the X-ray spectrum presenting overionized (recombining) 
plasma features. 

In recombining SNRs, the ionization temperature is 
higher than the electron temperature. This overionized 
plasma is a signature of rapid electron cooling. The ex¬ 
istence of recombining plasma (RP) was discovered in 
IC443 and W4 9B during the X-ray studi es of ASCA on 
six MM SNRs ([Kawasaki et al.|2002l|2005|). Additionally, 


the X -ray Imaging Spectrom eter, XIS, flKo yama et al. 
2007| onboard Suzaku (|Mitsuda et al.||2007|) has discov¬ 
ered strong radiative recombination continuum (RRC) 
features from the following MM SNRs: IC443, W49B, 
G35 9.1-Q.5, W28, W44 G346.6-Q.2, 3C 391, and G290.1- 
0.8 ([Y amaguchi et al. 2009 Ozawa et al. 20091 Ohnishi 


et al. | ^011 H bawada fc koyaiba||^ U\ 1 Uchida etHf2 012; 

Yamauchi et al.||2013| lErgin et al.||2014| IKamitsukasa etl 


there are other clues which may support the IVHVl-nature 


al.|[2014)_ 

Recently, |Yasumi et al.|(|2014|) studied G349.7+0.2 us¬ 
ing Suzaku data. 'They found that the spectrum of the 
whole SNR was well described by two optically thin 
thermal plasmas: a low-temperature (^0.6 keV) plasma 
in collisional ionization equilibrium (CIE) and a high- 
temperature (^1.24 keV) plasma in NEI. The low tem¬ 
perature plasma has solar metal abundances suggesting 
that the emission is dominated by the interstellar mate¬ 
rial. The high temperature plasma has super-solar abun¬ 
dances implying that the emission is influenced by the 
shocked ejecta. By using t he ab undance pattern of the 
ejecta component |Yasumi et al. (2014) estimated that 
the plasma originates from the core-collapse supernova 
(SN) explosion with the progenitor mass of ^35—40 M 0 
for this remnant. 

In this paper, we investigate the characteristics of the 
continuum radiation; thermal bremsstrahlung continuum 
and RRC in two regions, the brighter eastern (E) region 
and the western (W) region, by utilizing the superior 
spectral capabilities for diffuse sources of XIS onboard 
Suzaku. We report on our analysis of the E and W re¬ 
gions of G349.7+0.2. We analyze the GeV gamma-ray 
data, and probe the source morphology and variability 
of G349.7+0.2. In addition, we perform a detailed mod¬ 
eling of the GeV gamma-ray spectrum to understand if 
the gamma-ray emission from G349.7+0.2 might be the 
result of hadronic interactions between the SNR shock 
and the associated MC. 

2. DATA ANALYSIS AND RESULTS 
2.1. X-Rays 

2.1.1. Observation and Data Reduction 

G349.7+0.2 was observed with Suzaku on 2011 Septem¬ 
ber 19 with XIS for a total exposure time of 160 ks (obser¬ 
vation ID: 506064010). The XIS consists of three active 
X-ray CCD detectors (XIS 0, 1, and 3). XIS0 and XIS3 
have front-illuminated (FI) CCDs, while XIS1 contains 
a back-illuminated (BI) CCD. 

For data reduction, we used the High Energy Astro¬ 
physics Software (HE A Soft) package version 6.16 and 
the calibration database (CALDB) ve rsion 20140520. We 
used xspec version 12.8.2 (Arnaud 1996) for spectral 
analysis. The redistribution matrix an cl the auxiliary 
response functions were generated by xisrmfgen and 
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Figure 1 . Suzaku combined images of XISO, XIS1, and XIS3 
in the 0.3—10.0 keV energy band. The boxes are the spectral 
extraction regions correspo nding to the E and W of the SNR. 
CXOUJ171801.0—372617 (Lazendic et al. 2005) is indicated 
as a black cross. 


xissimarfgen ( Ishisaki et al.| 2007), respectively. 


2.1.2. Estimation of background 

Since G349.7+0.2 is located on the Galactic plane, the 
background of this SNR consists of the instrumental non- 
X-ray background (NXB), the Cosmic X-ray background 
(CXB), and the Galactic ridge X-ray emission (GRXE). 
Considering the position dependence of the GRXE, we 
constructed the background spectrum from the nearby 
blank sky (observation ID: 1000280 30) as the backg round 
is dominated by the GRXE (|Uchiyama et al.|2013|) . The 
NXB spectra for the source and background were esti¬ 
mated from the Suzaku database of^dark-Earth obser¬ 
vations using the procedure of Tawa et al. (2008|. The 
source and background spectra were made by subtract¬ 
ing the NXB spectra using mathpha. NXB-subtracted 
background spectrum was subsequently subtracted from 
the source spectrum using XSPEC. 

2.1.3. Spectral analysis 

Figure [l] shows Suzaku combined XISO, XIS1, and 
XIS3 image of G349.7+0.2 in the 0.3—10.0 keV en¬ 
ergy band. In order to study the spectral properties, 
we extracted the data from a 1'7 x 1.'3 size box cen¬ 
tered at R.A.(2000) = 17 h 18 m 03!7 and decl.(2000) = 
—37° 26'56'/2 for the E region and from a same-size box 
centered at R.A.(2000) = 17 h 17 m 5778 and decl.(2000) = 
—37°25 / 37 / /9 for the W region. These regions are shown 
with black rectangles in Figure [l] The E region corre- 
sponds to regions 1, 2, and 3 in Figure [T| of Eazen d ic et| 
al. (2005) and the W region includes the regions 4, 5, 
and 6 in the same figure. 

For the E region, we first tried to fit the spectrum 
to a single-temperature plane-parallel shock model with 


0.1 



2 5 

Enorgy (keV) 


Figure 2. Suzaku FI spectra shown in the 1.0—8.0 keV en¬ 
ergy range with the absorbed V VP SHOCK plus RRC model 
and the gaussian line for the E region (upper panel) and the 
VVPSHOCK model for the W region (bottom panel). 


variable abundances for all elements (Z < 30) (VVP- 
SHOC K in xspec) modified by TBABS absorption 
model ([Wilms e t al.||2000|) in xsp ec. The abundances 
are set to wilm ([Wilms et al. 2000). In this fitting, the 
absorbing column density (Ah), upper limit of the ion¬ 
ization timescale (r u ), electron temperature (fcT e ), nor¬ 
malization, and abundances of Si, Ca, and Fe were left 
as free parameter s. The o the r elem ents were fixed at 
their solar values ( Wilms et al.p OOO). Although this fit 
gave an acceptable x 2 /degrees of freedom (dof) value of 
448/429, residuals at ^3.3 keV and ^4.1—4.5 keV re¬ 
mained. The first residual around 3.3 keV corresponds 
to Ar Ly<a line. The other residual could be caused by 
the RRC corresponding to the H-like Ar. 

Thus we added a Gaussian and RRC (REDGE in 
XSPEC) components to the model. The line energy of 
the Ar Lya was fixed at 3.3 keV, while the normaliza¬ 
tion was a free parameter. We note that the Gaussian 
line width parameter was fixed to zero eV. For the RRC 
component, the energy was fixed at 4.4 keV. After adding 
the Gaussian and RRC components to the model, the re¬ 
duced x 2 value improved to 1.01 for 427 dof, which yields 
an F-test probability of 4.2 x 10 -4 . 
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Table 1 

Best-fit parameters with corresponding errors at 90% confidence level for the E and W spectra of G349.7+0.2. The spectral 

fits were carried out in the 1.0—8.0 keV energy band. 


Component 

Parameters 

E 

W 

TBABS 

VVPSHOCK 

N h (10 22 cm" 2 ) 
kT e (keV) 

Si (solar) 

Ar (solar) 

Ca (solar) 

Fe (solar) 
r u (10 11 cm -3 s) 

Norm (10 2 photon cm -2 s 1 ) 

q 7+0.3 

0.5 

1.2+0.1 
i 4+ 0 - 3 

1 (fixed) 

2 0 +0 ' 7 
Z > U -0.5 

o q+2.1 

*"*- 1.5 

13.916.1 

5.2ir 7 

9.5+0.3 

1.3+0.1 

1.2+0.1 

1.2+0.2 
i q +0.4 

1.5+0.4 
4.8i- 
5.8+0.3 

Ar Lya 

E (keV) 
o- (keV) 

Norm (10 —5 photon cm -2 s —x ) 

3.3 (fixed) 
0.0 (fixed) 

1 4+0-7 
i - 4 -0.8 


RRC H-like Ar 

E (keV) 

Norm (10 —5 photon cm -2 s —x ) 

4.4 (fixed) 
3 -°1 i'.5 



X 2 /dof 

432/427 

539.7/529 


reduced y 2 

1.01 

1.02 


We also applied RP model VVRNEI in xspec for E 
region. VVRNEI is a non-equilibrium recombining col- 
lisional plasma model with variable abundances for all 
elements (Z < 30). It is characterized by a constant 
electron temperature ( kT e ), initial temperature (&Tmi t ), 
and a single ionization parameter. Nh, kT e , the abun¬ 
dances of Si, Ca, Fe, r, and the normalization were free 
parameters during the fitting. We tried fitting by fixing 
/cTinit to a range of values (2.0, 5.0, 10.0, and 50.0 keV) 
and for each fit we obtained similar x 2 /dof and r values 
for each of the kT- in i t values. This model gave a x 2 /dof 
of 520.6/431, which required a large ionization timescale 
(r ~ 2.8 x 10 13 cm -3 s) showing that G349.7+0.2 is in 
ionization equilibrium and not in RP phase. 

As a next step, we applied an absorbed VVPSHOCK 
model for the W region. During the fitting, the absorp¬ 
tion, electron temperature, ionization timescale, abun¬ 
dances of Si, Ar, Ca, and Fe, and normalization param¬ 
eters were allowed to vary freely, whikythe other ele¬ 
ments were fixed at their solar values ([Wilms et al. 2000). 
We obtained an acceptable fit with a x 2 /dof value of 
539.7/529 and found no significant residuals. So, the 
spectral analysis of the E and W regions showed that 
there is no overionized plasma in G349.7+0.2. The spec¬ 
tral fits for the E and W regions are given in Figure [2] 
and the best-fit parameters are summarized in Table [if 
where all the errors are given at 90% confidence level. 

2.2. Gamma Rays 

In this work, Fermi -LAT data from 2008 August 
4 to 2014 March 27 were analyzed. The events- 
data were selected from a circular region of interest 
(ROI) with a radius of 18° centered at the SNR po¬ 
sition of R.A. (J2000) = 17 h 17 m 59M, deck (J2000) = 
—37°25 / 59 / /88 (R.A. = 259°.496 and deck = -37°.433). 
Using gtselect of Fermi Science Tools (FST), we chose 
the Fermi -LAT Pass 7 events suggested for the galac¬ 
tic point source analysis. We only selected the events 
coming from zenith angles smaller than 105°, to reduce 
the contributions from the albedo gamma rays from the 
Earth’s limb. 

To study G349.7+0.2’s spatial and spectral character¬ 
istics, the maximum likelihood fitting method (Mattox 



Figure 3. Gamma-ray TS map of G349.7+0.2 and its neigh¬ 
borhood with a bin size of 0°.01 x 0°.01. The black cross and 
circle are representing the GeV source location and its posi¬ 
tional statistical error from the 2nd Fermi -LAT catalog cor¬ 
responding to G349.7+0.2, respectively. The red cross and 
circle is the best-fit location and its statistical error circle 
(95% confidence level) found in this analysis, respectively. 
The blue contours show the Suzaku data from Figure [T] where 
the 3 contours represent X-ray counts higher than 790, 365, 
and 540 counts (cts). 


et al 
trally 


[T996) was employed on the spatially and spec- 
binned data within a square region of ~20° x 20° 
using the P7SOURCE_V6 version of the instrument re- 
spon se function. The analysis was p erformed using point¬ 
like (Kerr 2011 Lande et al. 2012|) (FST-v9r32pO) and 
the standard binned likelihood analysis (FST-v9r2 7pl) 
packages, both based on gtlike (Abdo et al. [2009|), for 
cross-checking the validity of the results. 

The model of the analysis region contains the dif¬ 
fuse background sources and all the point-like sources 
from the 2nd Ferrm-LAT source catalog located within 
a distance of 18° from the ROI center. We fixed 
all parameters of the point-like sources in the model, 
except 2 sources within the distance of 2° from the 
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ROI center, where we set the normalization and spec¬ 
tral parameters free. The standard diffuse background 
model has two components: the diffuse Galactic emis¬ 
sion (gal-2yearp7v6-v0.fits) and the isotropic compo¬ 
nent ( iso-p7v6source.txt ), which is a sum of the extra- 
galactic background, unresolved sources, and instrumen¬ 
tal background. The normalization of the isotropic com¬ 
ponent was kept fixed, while the normalization of the 
galactic diffuse background was set free during the anal¬ 
ysis. 

The test statistics (TS) parameter having larger values 
indicates that the null hypothesis (maximum likelihood 
value for a model without an additional source) is incor¬ 
rect, where the detection significance of a source can be 
approximated to the square root of the TS value. 


2.2.1. Location and Spectrum 

We searched for the best-fit location of G349.7+0.2 
within the ROI, which was found as R.A. (J2000) = 
259°.509 ± 0°.013 and deck (J2000) = -37°.455 ± 0°.013. 
The TS value of the best-fit positiorQ was enhanced by 
2.7 a over the position of 2FGL J1718.1—3725 in the 2nd 
Fermi -LAT catalog. Then the model was refitted using 
the b est-fit position to compute the TS map (Section 
2.2.2) and the spectrum. 

To check the functional form of the spectral energy dis¬ 
tribution (SED), we considered G349.7+0.2 as a point- 
like source. First, the power-law (PL) function was fit 
to the data between 200 MeV and 300 GeV. Because the 
spectrum deviates from a PL function, we tested if the 
gamma-ray emission is better described by a log-parabola 
(LP) or a broken power-law (BPL) function, where their 
functional forms are as follows: 


• Log-parabola: 

F(E) ip = N 0 (E/Eft) ( r i + r 2 ln (E/Ei,)) 

• Broken Power-law: 

F(E) bpl = N 0 (E/E 6 )" ri for E < E b 
= N 0 (E/E 6 )- r2 for E > E b 

Here, F(E) LP and F(E) PPP are representing the flux 
for the LP and BPL type spectrum, respectively. The 
spectrum breaks or turns at a certain energy of E&, where 
the spectral index changes from Ti to IV N 0 is the 
normalization parameter. 

Using different functions in fitting the spectrum of 
G349.7+0.2, the likelihood ratio, TS, was used as a mea¬ 
sure of the improvement of the likelihood fit with respect 
to the simple PL. The results to the spectral fits are sum¬ 
marized in Table [21 where the highest TS value (168) was 
found for the BPL fit. 

The PL resulted in spectral index of Ti = 2.00 ± 0.03, 
which is in agreement with the best-fit power-law index 
value given for 2FGL J1718.1—372 5 in th e 2nd Fermi- 
LAT catalog (~1.98), |Nolan et al.| fl20 12|. This result 
also m atches to the results obtained by | Castro fc Slane| 
(2010). The best-fit parameters for the BPL fit are N 0 

-li auw— i nwi—2 o-i 


= (2.30 ± 0.21) x HT 11 MeV -1 cm“ 2 s“\ Ti - 4.98, 
and T 2 = 2.28 db 0.04, where the given uncertainties are 
statistical. The total energy flux was found as (3.70 =b 
0.10) x 10 -11 erg cm -2 s _1 . 


1 R.A. (J2000) = 17 h 18 m 02M6, decl. (J2000) = -37°27T8700 


Table 2 

Spectral fit parameters for PL, LP, and BPL between 200 
MeV and 300 GeV assuming G349.7+0.2 as a point-like 


Spectral 

Model 

Photon Flux 
[ 10~ 8 ph cm -2 s -1 ] 

r 1 

r 2 

E h 

[MeV] 

TS 

PL 

3.59 ± 0.21 

2.00 ± 0.03 

- 


137 

LP 

1.66 ± 0.13 

1.67 

0.19 ± 0.31 

1410 

151 

BPL 

1.24 ± 0.14 

4.98 

2.28 ± 0.04 

568 

168 



Figure 4. The gamma-ray spectrum of G349.7+0.2. The 
Ferrm-LAT spectral data points are represented by the red 
filled circles with their corresponding statistical and sys¬ 
tematic errors shown in black and red color, respectively. 
The models of emission are shown by the thick blue (71°- 
decay component) and the dashed-dotted magenta (the 
bremsstrahlung emission component) lines. The parameters 
used to estimate the emission spectra are explained in Section 


We also calculated the systematic errors that stem 
from the uncertainties in the Galactic diffuse background 
intensity, where we foll owed the methods d escribed in 


Abdo et al. (2009a) and Castro et al. (2013). We varied 
the normalization value of the Galactic background by 
±6% from the best-fit value and used these new frozen 
values of the normalization parameter to recalculate the 
SED of G349.7+0.2. The systematic errors on the SED 
are shown in Figure [4] with red solid lines on top of the 
statistical errors given in black lines. 


2.2.2. Extension 

To investigate the morphology of G349.7+0.2, we cre¬ 
ated a 2° x 2° TS map of G349.7+0.2 and its neigh¬ 
borhood with a bin size of 0°.01 x 0°.01. The TS map 
shown in Figure [3] was produced with pointlike using a 
background model file, which contained all the point-like 
sources and diffuse sources, but excluded G349.7+0.2 
from the model. So, it shows the TS distribution of 
gamma rays originating dominantly from G349.7+0.2. 
In Figure [3] the blue contours represent the Suzaku com¬ 
bined XIS image in the 0.3—10.0 keV energy band (from 
Figure [l]) and the black cross and circle represent the lo¬ 
cation and its statistical error of the GeV source from the 
2nd Fermi -LAT catalog corresponding to G349.7+0.2, 
respectively. The best-fit location of the gamma-ray 
emission is slightly offset from the center of the X-ray 
remnant, as shown in Figure [3] by the red cross and cir- 
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cle, and it is closer to the south-east region of the SNR. 

To search for the energy dependent morphology, we 
split the data set into two energy ranges (200 MeV—1 
GeV and 1—300 GeV) and analyzed the data in each 
range. We found no significant gamma-ray excess at the 
location of G349.7+0.2 for the energy range between 200 
MeV and 1 GeV, but G349.7+0.2 was detected in the 
higher-energy range of 1—300 GeV with a significance of 
~lla using the BPL spectral model. 

Assuming that G349.7+0.2 has a PL/BPL type spec¬ 
trum, we checked the extension of G349.7+0.2 using the 
TS of the extension (TS ext ) parameter. TS ext is the 
the likelihood ratio comparing the likelihood for being 
a point-like source to a likelihood for an existing exten¬ 
sion. The pointlike code calculates the TS ext by fitting 
a source first with a disk template and then as a point- 
like source. In pointlike , the extended source detection 
threshold is defined as the source flux, where the TS ext 
value averaged ove r many statistical realizations is 16 
(Lande et al. 2012). Simulation studies showed that to 
resolve a disk-like extension with a radius (r) of 0°.l at 
the detection threshold, the source with a spectral index 
value of 2.0 and 2.5 must have a minimum flux of 3 x 
10 7 p h cm 

tively (Lande et al. 2012). 


s- 1 and 2 x 10“ 6 


ph cm 2 s x , respec- 


Table 3 

Fit results of disk-like extension model applied to G349.7+0.2 
gamma-ray data. 


Spectral 

Model 

Longitude 

[°] 

Latitude 

[°] 

Sigma 

[°] 

PL 

BPL 

349.70 ± 0.01 
349.70 ± 0.01 

0.158 ± 0.013 
0.154 ± 0.013 

0.003 

0.005 


For G349.7+0.2 having either a PL or a BPL type spec¬ 
trum, the TS ext value was found to be about zero, after a 
disk template fitting, where the fit parameters are shown 
in Table |U This indicates that a disk-like extension with 
r ^ 0°. 1 could not be resolved at the integrated flux level 
of 3.59 x 10 -8 ph cm -2 s -1 for a PL type spectrum with 
a spectral index of 2.0 and of 1.24 x 10 -8 ph cm -2 s -1 
for the BPL type spectrum with a spectral index of 2.3 
(T 2 ), respectively. 

2.2.3. Variability and Pulsation 

We first looked for long term variability in the light 
curve of G349.7+0.2 by taking data from the circular re¬ 
gion of 1° around the best-fit position. Figure [5] shows 
the 1-month binned light curve obtained after applying 
Fermi- LAT aperture photometry , where we checked for 
possible variations in the flux levels. If any or some of 
the flux data points are above 3cr, this would be an indi¬ 
cation of a significant variability (e.g. flare) in the circu¬ 
lar region of interest. In Figure [5] most of the flux data 
points remain within the la and 3a bands. Although 
one flux data point is above 3a, we note that it’s error 
bars are large, showing a large statistical uncertainty in 
this time interval. Thus, we conclude that there is no 
long term variability observed in the close neighborhood 
of G349.7+0.2. 

We have also checked if the spectral shape of 
G349.7+0.2 fits to the standard spectrum of a pulsar, 
Power Law with Exponential Cutoff (PLEC). The best- 


fit cutoff energy was found as 24.60 ± 4.95 GeV, which 
is an order of magnitude away from the r ange of typi¬ 
cal pulsar cutoff energies (Abd o et al.|201 0). The PLEC 
fit did not show a significant improvement over the PL, 
BPL, and LP spectral fits. 


6.50 

6.00 


w 5.50 

'§ 

1 5.00 
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Figure 5. The 1-month binned light curve showing no signif¬ 
icant variability in G349.7+0.2. 
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2.2.4. Molecular Environment 

We used th e CO d ata ta ken by the 12 m telescope of 
the NRAO (|Reynoso & Magnum 20011 with a good an¬ 
gular resolution (54") to understand the molecular en¬ 
vironment around G349.7+0.2. We integrated the CO 
spectrum over the whole range of velocities from —137 
to +57 km s -1 to obtain the velocity integrated CO in¬ 
tensity (Wco)- The W co averaged over the region with 
a radius of 0°.025 covering the X-ray remnant was found 
to be ^150 K km s -1 . Using the CO-to-H 2 c onversio n 
factor of X = 1.8 x 10 20 cm -2 K -1 km s -1 (Dame et 
al. 200l|), we found N(H 2 ) = 2.7 x 10 22 cm . Using 
the IN (H i) value found by|Tian & Leahy (2014), we cal¬ 
culated the total column density as JN(total) = N(Hl 
+ 2N(H 2 ) = 7.12 x 10"“ 


22 cm 2 


Lazendic et al. 

2*2 


(2005) 


cm 


^ 22 cm 2 as- 


calculated the total column density as 7.3 x 10 2 
adding the contributions from Hi, 2.8 x 10 2 
suming T s = 140 K, and H 2 , 4.6 x 10 22 cm -2 . 

Figure [6] shows the CO map produced in the velocity 
range of {+14, +20] km s -1 , where white contours rep¬ 
resenting the TS values shown in Figure [3j Also shown 
with black plu s m arkers are the five masers measured 
by Frail et al. (1996), which are located inside the blue 
Suzaku X-ray contours of 540 cts. 

The progenitor of the G349.7+0.2 has exploded within 
a large ringlike CO shell (seen as the color-coded arc 
structure in Figure |6|, a remnant of an older SN explo¬ 
sion containing denser clumps of MCs. The exp losio n 
probably happene d within the inter-clump gas (Reynoso 
& Magnum 2001) of the CO shell, at a location between 
the edge of the CO shell and an MC clump inside this 
shell. In our interpretation we assume that while the 
SNR’s shell was expanding, a small section of it broke 
out of the CO shell boundary and entered the rarefied 
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Figure 6. The CO data (Reynoso V Magnum 2001) inte¬ 
grated in the +14 to +20 km s -1 velocity range around the 
G349.7+0.2 location, where the black contours are for W co 
at 10, 18, 24, and 31 K km s _1 . The white contours represent 
the gamma-ray TS values of 90, 110, and 130. The blue con¬ 
tours show the X-ray counts from Suzaku , which are at 190, 
365, and 540 cts . The black plus markers are the five masers 
found by Frail et al. (1996). 


interstellar medium (ISM) moving in a direction toward 
the observer. 

The X-ray morphology of G349.7+0.2 is similar to the 
radio image consisting of a two-ring structure, where one 
ring is smaller and brighter than the other. The smaller 
ring might be the part of the SNR shock front enter¬ 
ing into a higher density region , an MC clump called as 
‘cloudl’ by Reynoso & Magnum (2001), and interacting 
with it. The shock of G349.7+0.2 entering ‘cloudl’ was 
confirmed to be of type C, wh ich allows OH 1720 MHz 
masers to form (Lockett et al.| 1999). The larger ring 
structure represents the S1NR shocFexpanding into the 
relatively less denser medium of the CO shell. 

Therefore, we calculated the proton density for two 
regions: one having a radius of 0°.015 corresponding to 
the inner smaller ring, and another one with a radius of 
0°.022 representing the larger outer ring. The smaller 
ring radius was chosen such that it covers the highest 
number of X-ray counts (shown with the innermost blue 
contour for >540 cts in Figure [6|. This region coincides 
mostly with the E region of the SNR selected in the X-ray 
analysis (Section 2.1.3). The larger ring was chosen to 
include the X-ray contour representing the X-ray counts 
higher than 365, where the W region is also partially 
covered. Both rings are within the lOcr contours of GeV 
gamma-ray TS distribution. To calculate the average 
proton density within the regions surrounded by these 
two rings, we first estimated the size of the emission re¬ 
gions to be ^3 and ^4.5 pc using the distance to the SNR 
to be 11.5 kpc. Then assuming a spherical geometry of 
the emission region, the proton density (n) correspond¬ 
ing to the inner and outer ring region was found to be 
81 and 35 cm -3 , respectively. 


2.2.5. Modeling and Interpretation 

To understand whether the leptonic or hadronic sce¬ 
nario is dominating the SED of G349.7+0.2, we started 
the multi-wavelength modeling assuming a BPL type of 
spectra for both the electron and proton distributions 
inside the source: 


dN t 


e,p 


dE 


N e,p E e,p for E e,p < E br 


E P 


E r . 


= Nl p E-$ exp I -+^ ) for E br < E ej> < E c (1) 


where E e , p is the electron or proton energy and E^ r is the 
spectral break energy with the spectral index changing 
from a to (3. E c is the maximum energy of electrons or 
protons. N \ p and Ng are normalization constants. 

For the multi-wavelength modeling we used the radio 


data 

(Shaver et al.|1985 

Whiteoak & Green 11996 Green 

2009 

and the gamma-i 
5.S. The H.E.S.S. fl 

:ay data from Fermi -LAT and 

ux data points (Abramowski et 
rgenta color in Figure [7| During 

al. 2014) are shown in m; 


and 2 TeV for electrons. The value of E c was chosen such 
that it would not cause the spectra of bremsstrahlung 
and inverse-Compton (IC) overestimate the TeV fluxes, 
when they are considered the dominant models for ex¬ 
plaining the data. Similarly, E c for protons was chosen 
such that the hadronic model does not overestimate the 
TeV fluxes. 

We first considered a purely leptonic model to explain 
the observed fluxes from radio wavelengths to gamma 
rays. In the case of lep tonic m odel, bremsstr ahlung and 
IC emission processes (Blumenthal & Go uld||19 70) were 
used for explaining the GeV—TeV dataT The total ener¬ 
getics of the electrons and protons need to be known to 
estimate the electron to proton ratio, defined as the ra¬ 
tio of dNe/dE and dNp/dE at E = 10 GeV, and the 
different model contributions. We first estimated the 
bremsstrahlung spectrum using a fixed value of ambient 
proton density. Parameters of the leptonic model then 
were adjusted in such a way that the bremsstrahlung 
spectrum could explain the observed fluxes at GeV—TeV 
energies. These parameters then were fixed and they 
were used in computing the IC and synchrotron spec¬ 
tra of the leptonic model. To explain the observed radio 
fluxes, the magnetic field was adjusted to get radio syn¬ 
chrotron spectrum. Using the same ambient proton den¬ 
sity we estimated the hadronic contribution to gamma- 
ray spec trum resulting from the decay of i r° (Kelner et 
aTp0 06). Secondly, we repeated a similar procedure for 


tEe case of IC process as the dominant model explaining 
the observed fluxes at GeV—TeV energies. And at last, 
we considered the 7r°-decay model dominating among ajl 
the pr ocess es, which is a purely hadronic model (Kelner 
et al. 2006). The parameters used for all three dominaT 


mg processes are shown in Table [4] for the two ambient 
proton densities of 35 and 81 cm -3 derived in Section 

f2~2~n 

Since the ambient proton density of 81 cm -3 repre¬ 
sents the interaction region corresponding to the smaller 
and brighter SNR shell, including the five OH masers, 
the modeling results for this density will be important 
in terms of giving us clues about the proton spectrum. 
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Table 4 

Parameters of multi-wavelength model. For bremsstrahlung and 
7r°-decay dominated models ambient proton density is considered 
to be 35 cm -3 and 81 cm -3 . 




The best-fit parameters for the proton spectrum were ob¬ 
tained by a y 2 -fitting procedure to the flux points. The 
dominating models of bremsstrahlung, IC, and 7r°-decay 
for the ambient proton density of 81 cm -3 are shown in 
the top, middle, and bottom panels of Figure [7| respec¬ 
tively. The estimated parameters are a = 1.89 db 0.14, 
P = 2.42 ± 0.03, and E hr = 12 GeV. The y 2 /dof is 
estimated to be ~ 0.53. The best-fit gamma-ray spec¬ 
trum resulting from the decay of i r° is shown in the 
bottom panel of Figure [7] with the red solid line. The 
estimated total energy can be written as W p ~ 3.13 x 
10 49 (81 cm _3 /n) erg, where n is the effective gas num¬ 
ber density for p-p collision. 


3. DISCUSSION 

G349.7+0.2 is an example of an ME SNR having an 
irregular she ll that lacks c entral cavity in radio wave¬ 
lengths (Shaver et al. 1985]), where the radio morphology 
is very similar to the one in X-rays. In both wavelengths 
the E side of the SNR is brighter. The most likely expla¬ 
nation for G349.7+0.2 not exhibiting a typical MM mor¬ 
phology is because we are likely seeing this SNR expand¬ 
ing into a density me dium with a large dens ity gradient 


of an angle of ~45° (|Lazendic et al. 2005). The pro¬ 


genitor (~35—40 M 0 ) of G349.7+0.2 possibly exploded 
within the ring-like CO shell having an average density 
of 35 cm -3 , where the location of the explosion may be 
between a dense MC clump and the boundary separat¬ 
ing the CO shell from the ISM. The radio continuum 
morphology shows two rings, where the smaller ring rep¬ 
resenting the part of the shell expanding into the higher 
density medium of the MC clump. While this part of 
the SNR shell might be expanding in the direction away 
from the observer, it is interacting with the MC clump 
producing 1720 MHz OH maser emission. The side of 
the SNR shell expanding towards the observer, on the 
other hand, is moving into a less dense region inside the 
CO shell and a smaller section of it might break out of 
the CO shell into the ISM. 

The stellar wind of the progenitor forms a CIE plasma 
in the dense circumstellar medium (CSM). When the 
blast wave breaks out of the dense CSM and expands adi- 
abatically into a cavity or into the rarefied I SM, it causes 
rapid cooling and the formation of the RP (|Itoh fc Ma- 


sai|1989||Moriya|2012| |Shimizu, Masai, & Koyama 2012|). 
So, we would expect the electron cooling for G349.7+0.2 
to happen through the adiabatic expansion mechanism, 
when a part of the shell closer to the observer has ex¬ 
panded into a lower density CO shell and then broke 




Figure 7. The different emission models fit on the multi¬ 
waveband data of G349.7+ 0.2 for the proton density yalue of 


81 cm 3 . The radio data (Whiteoak & Green 


1996 


Green 


2009) given in magenta filled triangles, are tit by the syn¬ 


chrotron emission model shown in dashed blue line. Fermi- 
LAT data points are represented by the blue filled circles and 
their corresponding s tatistical errors are also shown in blue. 
The H.E.S.S. results (Abramowski et al.||2014) are shown by 
the magenta filled squares. The solid light green and the 
double-dot-dashed magenta lines are the spectra correspond¬ 
ing to the bremsstrahlung and the IC models, while the solid 
red line shows the 7r°-decay contribution. The total emission 
is given in solid dark-green line. The top, middle, and bot¬ 
tom panel shows the bremsstrahlung, IC, and the 7r°-decay 
dominated model, respectively. 


out into the rarefied ISM. We searched for RRC in the 
Suzaku data from the E and W regions of G349.7+0.2. 
However, we did not detect any RP features in their spec¬ 
tra. The total column density calculated from the atomic 
and molecular gas data was obtained from a circular re¬ 
gion covering the X-ray remnant and it was found as 
N(total) = 7.1 x 10 22 cm -2 , where TVh value of the E 
and W regions (Table [l]) were found to be higher than 
the N(total) value. The column density value calculated 
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for the E region of G349.7+0.2 was found to be lower 
than the N h value for the W region, while the W region 
represents the portion of the SNR shell expanding inside 
the CO shell in other directions. The Suzaku analysis 
revealed the enhanced ejecta emission of Si, Ca, and Fe 
in the E region pointing to the reverse shock produced 
by the SNR-MC interactions. 

To understand whether the dominating gamma-ray 
emission model for G349.7+0.2 is leptonic or hadronic, 
we modeled all of the emission processes at GeV energies 
individually. We have adjusted the model parameters 
to explain the observed fluxes such that we get a set of 
parameters for each different dominant process. How¬ 
ever, some of the parameters estimated were either non¬ 
physical or inconsistent with the observed values. Based 
on these parameters, we could get some insight into the 
emission mechanisms inside G349.7+0.2. 

For the case of bremsstrahlung dominated model, the 
estimated electron to proton ratio was found to be much 
higher than the value (~0.01) observed on Earth. In the 
case of IC process, we found that the number of ambient 
proton density has to be less than or equal to 0.1, which 
is inconsistent with the measured proton density values 
reported in this paper. Moreover, electron to proton ratio 
was found to be much higher than the observed value. 
Therefore, both of these leptonic emission processes are 
insufficient to explain the observed gamma-ray fluxes. 

In the case of hadronic model, we considered observed 
electron to proton ratio to be about 0.01, and found that 
the hadronic model could explain the observed fluxes well 
without having any inconsistency in the model parame¬ 
ters. The magnetic field was found to be 150 /iG, which is 
close to the ma gnetic field v alue o f one of the OH masers 
measured by Brogan et al. (2000). 

The total gamma-ray luminosity was found as L = 1.58 
x 10 36 erg s -1 , similar to the first GeV-emitting SNRs 
that were d iscovered by Fermi- LAT, e.g . IC44 3 (pVbdo 



2009a), W44 (Abdo et 


and 3C 3yl(|Ergin 


interacting MM RNKs 

with gamma-ray luminosities higher than 10 35 erg s -1 . 
In order to explai n the reported GeV - TeV g amma-ray 


flux of G349.7+0.2, Abramowski et al. (2014) estimated 
that the total proton energy is unreasonab. 
an ambient proton density of 1.7 cm -3 


•ly large for 
As a result, 

they concluded that the TeV gamma-ray emission possi¬ 
bly cannot arise from the whole SNR shell, but from a 
region where the SNR is interacting with an MC. This 
conclusion also supports our work, which we have done 
in this paper considering a SNR-MC interaction. They 
also reported that a leptonic model is strongly disfavored 
when trying to explain the observed gamma-ray spec¬ 
trum. This results also agrees with our conclusion in 
this paper. 

One possible scenario for explaining the observed 
gamma rays is that G349.7+0.2 has evolved inside an MC 
having dense clumps. The SNR has primarily evolved in¬ 
side the inter-clump medium, compressing and shocking 
the dense material throughout its passage. During this 
process, it heats the material to X-ray temperatures, as 
well as accelerates the particles through diffuse shock ac¬ 
celeration. The particles can be either reaccelerated cos¬ 
mic rays trapped in the MC and the shell of the SNR or 


freshly accelerated protons entering the radiatively com¬ 
pressed MC. In this scenario, since the crushed clouds are 
thin, multi-GeV particles can escape from the shocked 
MC, which might be the reason of seeing a break in the 
proton sp ectrum (Bla ndford fc Cowie|1982 Uchiyama et 


ah]|2Q10 ). In an alternative scenario, escaping relativis- 


tic protons reach a nearby unshocked MC and produce 
7r°-decay gamma rays. For this scenario to happen, there 
must be GeV/TeV sources found outside the radio shell 
of G349.7+0.2_that_could produce these escaping pro¬ 


tons (Gabici et al. 


. However, there are no nearby 


cosmic ray sources to G349.7+0.2. 

All of the interacting SNRs having RP plasma, where 
most of them also presenting hadronic gamma-ray emis¬ 
sion, are of MM-type. Finding RRC features and gamma 
rays dominated by the hadronic emission would strongly 
suggest that G349.7+0.2 is an MM type SNR. We found 
that G349.7+0.2 has a gamma-ray spectrum, which is 
likely hadronic in nature, but lacks RP plasma. 


4. CONCLUSION 

G349.7+0.2 is expanding in a dense medium of MCs 
and emitting gamma rays by interacting with clumps 
of molecular material. We analyzed GeV gamma rays 
from G349.7+0.2 and found that the observed gamma- 
ray spectrum is likely hadronic in origin and is following 
the spectrum of parent protons, a BPL distribution with 
spectral index parameters of a = 1.89 and (3 = 2.42 and 
a spectral break at ^12 GeV. This suggests that protons 
are accelerated to high energies, possibly at the region 
of SNR shell interacting with the dense MC clump. Al¬ 
though the best-fit location of the gamma-ray emission 
is slightly offset from the center of the X-ray remnant, it 
is closer to the E region. Both the locations of the five 
OH masers and the whole E region of the X-ray remnant 
are within the 10cr contours of the GeV gamma-ray TS 
distribution. The low breaking energy of the proton spec¬ 
trum can be explained either by higher energy protons 
being already escaped through the crushed MC or the 
SNR shell at the earlier evolutionary stages of the SNR, 
or by the SNR shock expanding in a dense medium be¬ 
ing slowed down by this medium. We searched for RP in 
G349.7+0.2, whose MM-nature is still unclear, by study¬ 
ing the plasma structure of two of its regions using the 
Suzaku data. We found that the plasma in the E region 
is in ionization equilibrium, while the plasma in the W 
region is in the NEI state, showing that none of these 
regions are in the RP phase. 
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